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The relationship between vacancy ordering and magnetism in TlFei. oSe2 has been investigated via 
single crystal neutron diffraction, nuclear forward scattering, and transmission electron microscopy. 
The examination of chemically and structurally homogenous crystals allows the true ground state 
to be revealed, which is characterized by Fe moments lying in the afe-plane below 100 K. This is in 
sharp contrast to crystals containing regions of order and disorder, where a competition between 
c-axis and afc-plane orientations of the moments is observed. The properties of partially-disordered 
TlFei.eSe2 are therefore not associated with solely the ordered or disordered regions. This contrasts 
the viewpoint that phase separation results in independent physical properties in intercalated iron 
selenides, suggesting a coupling between ordered and disordered regions may play an important role 
in the superconducting analogues. 

PACS numbers: 74.70.Xa,75.25.-j 



The microscopic coexistence of superconductivity 
and magnetism at temperatures up to ~30K in 
Ai_j / Fe2- :c Se2 (A — K, Cs, Rb) has generated much 
interest, and a plethora of techniques have been em- 
ployed to understand the interesting properties of these 
materials. [H3] Structural studies on superconducting 
samples revealed antiferromagnetic order with a local 
magnetic moment of ^3/is/Fe, and a VEa x y/Ea super- 
structure of the ThCr2Si2 structure-type derived from Fe 
vacancy ordering. [51 [3] Other studies have focused on the 
electronic structure [SHI] and the superconducting pair- 
ing mechanism. [BHT3] as well as the influence of phase 
separation. [T3HT§] 

The coexistence of superconductivity and magnetism 
is now believed to be due to fine-scale phase separa- 
tion, with each property being associated with a different 
composition and/or degree of order. 16 21 Phase separa- 
tion over length scales of ~ 10- 100 nm has been observed 
via transmission electron microscopy (TEM) , [T3] HH ITS] 
nano- focused x-ray diffraction. |16j and is inferred from 
Mossbauer spectroscopy. [P7] Superconducting samples 
that are single phase (structurally and chemically ho- 
moegenous) have yet to be produced and thus the in- 
herent nature of the individual phases remains unclear. 
Flexibility in the chemical composition near the parent 
A0.sFex.6Se2 allows superconductivity to be realized. [22T 
[27] and may do so by promoting phase separation via lo- 
cal inhomogeneity. It has been suggested that the super- 
conducting phase does not have Fe vacancies. |19l |2"51 12"9"] 
Similar to other Fe-based superconductors, these super- 
conducting phases may be characterized by a doping level 
of approximately 0.15 electrons per Fe,[28l [30] such as 
Rbo.3Fe2Se2. The degree of order is also important. [19] 
131] however, and first principles calculations have shown 



that the vacancy and magnetic order greatly influence 
the Fermi surface. [5] 

Fe vacancies order into the ^/5a x y/ba supercell be- 
tween 460 and 580 K. At or slightly below this temper- 
ature, the spins align along the tetragonal c axis in a 
block-checkerboard antiferromagnetic structure (BCAF- 
c, where c indicates the alignment of the spins) . [H 0] The 
vacancy order is generally incomplete, however, with par- 
tial occupancies observed on at least one Fe site. [35] [33J 
Superconductivity, which is believed to exist in a phase 
that does not contain vacancy order, is only observed in 
samples that also possess regions of vacancy and mag- 
netic order. 

In contrast to the alkali-metal based compounds, 
TlFei.gSe 2 has a limited compositional window and the 
Tl site is always fully occupied, which reduces the de- 
gree of inherent disorder. While superconductivity is ob- 
served for mixed Tl / alkali-metal compounds, and was 
reported with a small volume fraction in TlFei 7Se9.|25j 
it has not been reproduced in TlFei 6 Sc 2 . This is per- 
haps associated with an improper electron count in 
TlFei. 6Se 2 . Based on the chemical formula, one would 
expect TlFei. 6Se2 to be a metal. However, TlFei. 6Se2 is 
observed to be an insulator. TlFei. gSe2 also displays dif- 
ferent magnetic behavior than the alkali-metal analogues. 

The BCAF-c structure is the only magnetic structure 
observed in the alkali-metal compounds. However, sin- 
gle crystals of 'partially-disordered TlFei. 6 Se 2 ' have been 
shown to possess additional magnetic phase transitions 
near 100 K and 140 K. [31] A detailed structural study re- 
vealed these transitions arc associated with a canting of 
the Fe moments toward the a6-plane, which only occurs 
between ~100 and 150 K. [35] These crystals are char- 
acterized by crystallographically coherent regions of or- 
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dered and disordered Fe vacancies. [35] In the absence of 
single-phase crystals, it is impossible to know if these 
magnetic transitions in partially-disordered TlFei.gSe2 
are inherent to the ordered or disordered regions, or are 
the result of some interaction between the two. 

In this Letter, we provide a detailed characterization of 
single crystal TIFei.eSe2 with complete chemical/ vacancy 
order. Such structural and chemical homogeneity has 
been elusive in ^4i_ !/ Fe2-a;Se2, and the presence of mul- 
tiple phases has complicated the structure-property in- 
vestigations in these complex materials. By obtaining 
ideal order, a previously unobserved spin reorientation is 
revealed with spins lying in the a6-plane for T < 100 K. 
This suggests a strong interaction between the ordered 
and disordered regions prevents this ground state from 
occurring in the partially-disordered crystals at low tem- 
peratures. 

Crystals with complete vacancy order were grown in 
melts of nominal composition TlFei.gSe2. In the same 
furnace, a nominal composition of TlFei.7Se2 produced 
crystals with magnetization behavior qualitatively sim- 
ilar to the 'partially disordered' crystals, which were 
grown from a melt of nominal composition TlFe 2 Se2.[31] 
Therefore, growth in iron-rich environments appears to 
inhibit the ordering of iron vacancies, perhaps due to mi- 
nor differences in Fc content that arc difficult to detect. 
In the fully ordered crystals, the vacancy ordering tem- 
perature is found to be ~463 K by differential scanning 
calorimetry (DSC), and an anomaly in the magnetic sus- 
ceptibility (x) was a l so observed at this temperature. An 
entropy change of — 7.4 J/mol-Fe/K was calculated from 
the DSC data, and this is ~75% of the entropy change 
expected for the simultaneous ordering of Fe vacancies 
and Fe magnetic moments. See the Supplemental Infor- 
mation] for additional details. 

The crystals were observed to be chemically homoge- 
nous and fully-ordered via aberration corrected scan- 
ning TEM (STEM) using high angular annular dark field 
(HAADF) imaging. The STEM image shown in Fig|jJ, 
was taken with the electron beam parallel to the [110] 
direction of the supercell. As a result of this orienta- 
tion, all of the Fe vacancies are aligned in columns par- 
allel to the electron beam and can be readily viewed as 
dark spots separated by four iron columns. A crystal 
structure drawing corresponding to the HAADF image 
is shown in Fig.[T]3 to highlight this ideal vacancy order- 
ing, and the atomic layers are labeled for ease of view- 
ing. Specimens viewed along [001] also revealed an or- 
dered and uniform vacancy superstructure, and a selected 
area diffraction pattern from this orientation is shown in 
Fig.[l]:. The selected area diffraction pattern shows two 
sets of V5a x y/Ea superstructure reflections, indicated 
by the vectors q and q', due to twinning along the (110) 
subcell planes. Images of the partially-disordered crys- 
tals can be found in Ref.[3SJ STEM data reveal the dis- 
ordered regions contain Fe vacancies (modulations in Fe 




FIG. 1: (color online) (a) Aberration corrected HAADF image 
of a chemically homogenous crystal, with electron beam along 
the [110] direction of the y/5a x y/Ea supercell, where all Fe 
vacancies are aligned in columns viewed as the dark spots. 
A corresponding crystal structure drawing is shown in (b). 
(c) Selected area diffraction with electron beam along [001] 
showing reflections from the y/5a x y5a superstructure with 
twinning yielding spots from two domains indicated by q and 
q'; the pattern is indexed according to the subcell. 



HAADF intensity), and thus these regions are not like 
the vacancy-free phases being reported in the supercon- 
ducting analogues. 

Refinements of the neutron diffraction data indicate 
complete occupancy of the Fel site (Wyckoff position 16i) 
and zero occupation of the Fe2 site (Wyckoff position 
4d), which yields the ideal v5a x V5a superstructure in 
TlFei.6Se2. In Fig.[2]3, the intensity of the vacancy order 
(020) peak is shown to be independent of temperature 
below ~460 K, indicating a saturation of vacancy order 
at all temperatures probed by neutron diffraction. 

The magnetic susceptibility of fully ordered TlFei.6Se2 
single crystals is shown in Fig. [2^,, where a phase transi- 
tion is clearly observed near 100 K. Above this transition, 
the susceptibility is largest when the magnetic field (H) 
is perpendicular to c, and for temperatures below the 
transition the susceptibility is largest for H || c. This 
indicates a switch in the easy magnetic direction, with 
moments parallel to c above — 100K and perpendicular 
to c below -100 K. 

The reorientation of the magnetic moment is con- 
firmed via single crystal neutron diffraction (Fig.[2]D,c) 
and synchrotron radiation based nuclear forward scat- 
tering (Fig. [3]). Figure [2}d shows the integrated intensity 
for three Bragg peaks observed via neutron diffraction 
on a —50 mg single crystal. The intensities of the two 
peaks with magnetic contributions change rapidly near 
100 K. The (121) peak intensity is associated with mo- 
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FIG. 2: (a) Magnetic susceptibility of oriented single crystals 
reveals a change in the easy magnetic axis near 100 K. (b) 
The integrated neutron diffraction intensity (in counts) for 
the vacancy-order (020) peak demonstrates a saturation of 
the vacancy ordering at high temperature. The integrated 
intensity of the (121) is proportional to the c-axis component 
of the moment, and the (105) peak gains intensity when the 
moment has an ab-plane component; together, these magnetic 
peaks reveal a change in magnetic orientation near 100 K. (c) 
Projections of the refined moments onto the a&-plane and c- 
axis for fully ordered crystals and in (d) the same is shown 
for partially-ordered crystals, (e) The c-lattice parameters 
obtained by x-ray diffraction. 



merits aligned with the c-axis and is beginning to saturate 
before it quickly diminishes near 100 K. At this temper- 
ature, the intensity of the (105) peak increases sharply 
due to the onset of the new magnetic order. The (105) 
peak gains intensity when the et&-plane component of the 
moment increases, as well as when the Fe layer develops a 
corrugation. |35j The moment reorientation is highlighted 
in Fig. [2]: where the refined moments are presented, and 
reach ~3/xs/Fe in both magnetic structures. A similar 
moment is obtained via first principles calculations for 
the BCAF-c (Supplemental Information), and the cur- 
rent calculations are consistent with previous reports. [B] 
The canting of the Fe moment in the partially- 
disordered TlFei.6Se2 between 100 and 150 K is shown in 
Figure [2ji435J A competition between the BCAF-c and 
an in-plane structure exists in the partially-disordered 
crystals, where the BCAF-c structure is present below 
100 K and above 150 K, and thus the corresponding mo- 
ment has no a6-plane component. When this canting is 
lost below ~100K, an abrupt increase in the c-lattice 



parameter is observed (see Fig. (2^).[3H 135] The mag- 
netic transition in the fully ordered crystals is not associ- 
ated with an anomaly in the c lattice parameter, though, 
and this is consistent with the smooth behavior in c near 
~T50K where the canting of the moment begins in the 
partially-disordered crystals. 

The moment reorientation is also observed via nuclear 
forward scattering (NFS),[3S] the time-analogue of Moss- 
bauer spectroscopy. NFS measurements were completed 
at beam line ID18 of the European Synchrotron Radia- 
tion Facility, [37] with the crystallographic c-axis parallel 
to the x-ray beam. The spectra were fit by a conventional 
routine. [35] 

The NFS data at temperatures slightly above and be- 
low magnetic transition at 100 K are shown in Fig. [3] 
These results clearly indicate the change of the hyper- 
fine structure across the transition, and confirm the mo- 
ment orientation changes from parallel to c above 100 K 
to perpendicular to c below 100 K. At 110 K, the data are 
well described by a model with AFM Fe moments aligned 
with the c-axis, yielding a magnetic field 14=26.1 T, a 
quadrupole splitting of AEq=1.18 mm/s, and the angle 
between main axes of those interactions is 46°. These 
values are in excellent agreement with those observed for 
Rb . 8 F ei . 6 Se 2 .[n] 

The NFS data at 90 K were fitted with 2 Fe sites with 
moments aligned in afe-plane, and the angle between them 
is utilized as a fit parameter. This model can generally 
describe any complex, non-collinear magnetic structure 
with moments in the plane. Independent of the model 
utilized, a change of the carrier frequency from that re- 
lated to the beats of I s * and 6^ lines of Mossbauer spec- 
trum (at 110 K) to that related to the beats of 2 nc ^ and 
5 lines (at 90 K) is observed, indicating a reorientation 
of the Fe moments. An appearance of the strong 2 nc ^ 
and 5^ lines indicates the a6-plane alignment of the Fe 
magnetic moments below 100 K. See the Supplemental 
Information for further details on NFS data collection 
and analysis. 

Below 100 K, the neutron diffraction data were mod- 
eled using several different magnetic structures. First, 
the eight possible magnetic structures within the ex- 
pected tetragonal symmetry were examined. Of these, 
the best fit (i?=0.052 at 5K) corresponds to the non- 
collinear magnetic structure (74/m) shown in Fig.[4^i 
with moments pointing 10.8° away from the vacancies. 
Fixing the moments towards the vacancies results in a 
similar quality fit. 

Alternative magnetic models require lowering the sym- 
metry of the magnetic structure and assuming equal 
twin fractions of collinear magnetic domains. The mon- 
oclinic model (72'/m) shown in Fig.|4]D (and a twin gen- 
erated by rotating 90 degrees about the c-axis) can de- 
scribe the data well {R = 0.045 with 2.99(4)^ B /Fe, at 
5K). Other possible twinned collinear magnetic struc- 
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FIG. 3: (color online) (color online) Time evolution of NFS 
for fully ordered TlFei.6Se2 crystal with the c-axis parallel 
to the x-ray beam measured above and below the magnetic 
transition. The data are fit by models (red lines) with mag- 
netic moments aligned with the c-axis at 110 K, and lying in 
the aft-plane at 90 K. The insets are the model simulations of 
Mossbauer spectra, with the lines numbered in the inset of 
(b). 




FIG. 4: (color online) The Fe (blue) lattice with differernt 
magnetic structures that describe the neutron diffraction data 
below 100 K. (a) A non-collinear magnetic structure that de- 
scribes the neutron diffraction data and retains the symmetry 
of the crystal structure, (b) An in-plane, block-checker anti- 
ferromagnetic structure that describes the neutron diffraction 
data provided an equal mixture of domains rotated 90° about 
c exists. In both (a) and (b), the moments (arrows) are anti- 
ferromagnetically aligned between Fe planes along the c-axis, 
and vacant Fe sites are shown as light grey. 

ture models were tested but no others gave reasonable 
agreement. 

Theoretical calculations have suggested non-collinear 
structures are stable in the vacancy-ordered systems.^? 
Collinear, in-plane structures are observed in the major- 
ity of Fe-based superconducting families, [39] but the 
collinear structure shown in Fig.[4j3 would be unique to 
this system. The measurements presented here are inca- 
pable of conclusively identifying the magnetic structure 
below 100 K. In the Fe-based superconductors, structural 
and magnetic transitions are typically coupled. There- 
fore, it seems reasonable to expect that if the magnetic 
symmetry were lowered to monoclinic, to allow the struc- 
ture in Fig.[4j}, then the nuclear structure would also 



distort. Such a distortion is not evident in the current 
diffraction data. As such, the preference would typically 
be for the higher symmetry structure (Fig. [4^,), which 
also provides a good description of the data. Regard- 
less of the exact spin structure, this reorientation of the 
spins near 100 K is clearly driven by the complete order- 
ing of vacancies, and the dominance of this ground state 
is unique to this phase-pure material. 

The transitions at 100 K in the fully-ordered and 
140 K in the partially-disordered crystals are similar. In 
both cases, the moments develop an a6-plane compo- 
nent upon cooling through the transition. The struc- 
tures in Fig. [4] were not observed via neutron diffrac- 
tion for the partially-disordered crystals, though the NFS 
data suggest the Fe moments lie within the et6-plane at 
110 K (see Supplemental Information). In either case, 
the BCAF-c order with spins along c is the ground state 
in the partially-disordered crystals below 100 K. There- 
fore, the disordered regions, which are likely to be non- 
magnetic, [17] have a strong influence on the observed 
magnetic behavior. Consistent with these results, it 
has been recently suggested that superconductivity is in- 
duced in KFe2Se2 when this KFe2Se2 stoichiometry in- 
teracts with the vacancy-ordered K0.8Fe1.eSe2 phase. [29] 

Given the magnetoelastic nature of these materials, [SI 
1331 |4"U] it is not surprising that disordered and ordered 
(non-magnetic and magnetic) regions would be strongly 
coupled. These regions are crystallographically coher- 
ent and approximately 10-20 nm in a6-plane and 2-4 unit 
cells along c,[35] which allows them to interact via lat- 
tice strain. A detailed analysis of the STEM data for 
partially-disordered TlFei. 6 Se 2 reveals the c lattice pa- 
rameter is approximately 2% larger in the disordered re- 
gions than in the ordered regions, and the change in a 
is less than ~0.3%. Therefore, modulations in c are the 
most likely source for the strain interaction. 

As demonstrated in this Letter, the absence of disor- 
dered regions allows TlFei. 6 Se 2 to complete a spin re- 
orientation at ^100 K, thereby unveiling the true mag- 
netic ground state associated with the y/Ea x \/5a su- 
percell. This demonstrates a strong influence of the 
secondary phase on the behavior of the primary phase, 
and thus an interaction between the ordered and disor- 
dered regions determines the properties of the partially- 
disordered crystals. These results have implications for 
the superconducting alkali-metal iron selenides, which 
also possess regions of order and disorder separated on a 
fine scale. Understanding the interaction between these 
regions may prove critical in identifying the true nature 
of superconductivity in those materials. 
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